Unusual oxygen isotope effects in cuprates — importance of doping 
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A recent angle resolved photoelectron spectroscopy ( ARPES) study by Douglas et alP on oxygen 
isotope exchanged Bi2Sr2CaCu2 0g+i superconductors reported an absence of isotope effect at op- 
timal doping, questioning the previous work by Here, we report a new result that sheds light 
on this puzzling discrepancy as well as the nature of the electron lattice interaction in the cuprates: 
the anomalous isotope effect at optimal doping^, re-confirmed here, vanishes on a mere 2 % over- 
doping of holes. This result implies a rapid change of the nature of the electron-lattice interaction 
near optimal doping. We also find that the data by Douglas et alP^ are actually characteristic of 
significantly over-doped samples, not of optimally doped samples as they claimed. 
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Whether there is a critical doping concentration at 
which the electronic structure of hole-doped high tem- 
perature superconductor material undergoes dramatic 
change has been a subject of debate over the past. For 
example, pump and probe experiments^ revealed a dra- 
matic change of the quasiparticle recombination dynam- 
ics near optimal doping. Other thermodynamic and 
transport data have also pointed toward the possibility of 
a critical doping concentration somewhere in the vicinity 
of optimal doping^'. In this paper we shed some light on 
this question by examining the effect of oxygen isotope 
substitution ( 16 0^ 18 0) on the photoemission data. Our 
main conclusion is that the isotope effect is very sensitive 
to even a minute amount of doping in excess of the opti- 
mal doping. This result suggests that the electron-lattice 
coupling undergoes a sharp change across the optimal 
doping. 

In this paper, we extend our previous study of the 
oxygen isotope effect of the angle resolved photoelectron 
spectroscopy (ARPES) on Bi 2 Sr 2 CaCu208+(5 to slightly 
over-doped samples. The new data on optimally doped 
samples show anomalous isotope effect in off-nodal re- 
gions, in good agreement with previous results^'. In 
addition, as x, the amount of doped holes per unit cell 
of the 2-dimensional Cu02 layer, increases by a mere 2 
%, the isotope effect disappears within the error bar of 
the experiment (w ±10 meV). These results are consis- 
tent with the abrupt c hange s observed by other probes 
at similar doping valuesSHEI. 

In light of these findings, we respond to a recent com- 
ment on our previous work by Douglas et alP, who re- 
ported the absence of an anomalous isotope effect for 
"optimally doped" samples. Indeed, judging from the 
published data of Douglas et alP we estimate the doping 
level of their sample to be at least 2 % in excess of op- 



timal doping, which explains why they do not observe a 
significant isotope effect. 

The new ARPES data reported here were obtained at 
beam line (BL) 7 of the Advanced Light Source (ALS) of 
the Lawrence Berkeley National Laboratory, using pho- 
tons of energy 100 eV at temperature w 25 K. The energy 
resolution was set to ~ 60 meV or « 35 meV while the 
angular resolution was set to w 0.35°. Note that although 
these energy resolution values are rather large compared 
to the best available energy resolution (1 ~ 10 meV), 
they are comparable to the resolution used in Ref. [TJ and 
are sufficient for discussing high energy peaks, the main 
focus of this paper, whose intrisinc widths are well above 
100 meV. Despite the relatively large energy resolution, 
the peak positions described below, extracted from peak 
maxima, have an uncertainty of only ±10 meV. We also 
discuss (in Fig.'s 3 and 4a) data obtained at BL 10 of 
the ALS, as described before^. It is important to note 
that both in the experiments rep orted here and in those 
reported in our previous worH 2 ^ the two oxgyen-isotope 
samples were measured one right after the other in an 
otherwise identical experimental set-up. 

Fig. [TJshows data obtained on optimally doped samples 
and over-doped samples, with x ~ 0.16 and x » 0.18 
respectively As we shall discuss later, although usually 
it is difficult to estimate the absolute doping level to an 
accuracy better than 1 % (i.e. 0.01 for x), ARPES is 
sensitive to a doping change of 1 %. 

Fig. [lja) shows a constant energy map at the chemi- 
cal potential (/i) for x = 0.16 sample, giving an overview 
of the Fermi surface. This map fits well with a tight 
binding parameterization^, whose Fermi surface and its 
replica arc plotted as solid lines. The high photon energy 
used here allows to enhance the super-structure replicas 
"SI" and "S2," first order and second order replicas due 
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FIG. 1: (a) Constant energy map at the chemical potential (fi) 
on optimally doped Bi2Sr2CaCu208+5 sample. The boundary 
for the first Brillouin zone is shown as a square, and four 
selected high symmetry points are explicitly indicated in unit 
of (tt / a,n / a), where a is the lattice constant of the CuC>2 
plane, (b) Angle integrated photoemission data obtained on 
three different sets of samples, at optimal doping (a; = 0.16) 
and a slight over doping [x = 0.18). Data are normalized to 
the area under the curve for E — n > —0.5 eV. For ease of 
view, each pair of curves has been shifted in intensity, which 
approaches zero at far right end of the energy axis for all 
curves. 



to the well-known surface reconstruction in the Bi-oxide 
layer. As we will discuss extensively below, these are 
key features for determining the doping with high ac- 
curacy (<C 0.01). It is important to note that, because 
of the geometry and photon energy used in this exper- 
iment, the ARPES signal is a local minimum along the 
nodal direction, consistent with previous findings using 
the same geometry^. This situation is very different from 
the previous one^ in which the nodal intensity was a local 
maximum. For this reason, in this paper we focus on the 
isotope effect only in the off-nodal region, where we have 
sufficient signal to noise ratios to discuss large isotope 
effectP. 

Fig. [TJb) shows angle integrated data as a function of 
the oxygen isotope for three pairs of isotope-exchanged 
samples and two different doping values. The upper two 
sets of data were obtained with a total energy resolution 
of ss 60 meV while the bottom set of data with a total 
energy resolution of ~ 35 meV. These curves present the 
main result of this paper: the isotope induced changes 
are strongly doping dependent. Note first that there is 
a clear isotope-induced shift of the broad high energy 
bump at —0.2 eV by « 30 meV at optimal doping (x = 
0.16), in agreement with previous reports^. Surprisingly 
the isotope shift is reduced to less than 10 meV when x 
increases by a mere Ax = 0.02. 

Douglas et al.^ have argued that the isotope-induced 
shift observed at optimal doping is an artifact caused by 
errors in sample alignment. We have chosen to compare 
angle integrated data^, as they are clearly immune to 
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FIG. 2: (a, b) EDCs for momentum value (0, 0.5), showing a 
finite isotope shift for optimally doped sample (a) and no shift 
for over-doped sample (b). (c) EDC peak positions compiled 
for three cuts, arising from main (cut 1) and superstructure 
replicas (cuts 2,3) for x = 0.16. When momentum values for 
cuts 2,3 are normalized, by subtracting relevant superstruc- 
ture wave vectors, it can be seen that all of these three cuts 
are defined on the same vertical line, passing through points 
(0, 0) and (0, 1), as shown in the inset. The horizontal axis of 
this plot is defined by this normalized momentum value. Note 
also that the graphs for cuts 2,3 are shifted in the increments 
of -0.2 eV for easy view of data, but it should be understood 
that each cut has an approximate energy range -0.2 eV to -0.4 
eV. (b) EDC peak positions for x = 0.18, plotted in the same 
manner as in (c). 

such systematic errors. Fig. [T] shows that the isotope 
induced shift is indeed an intrinsic effect and that the 
disparity between our previous data and those of Rcf. 1 
can not be attributed to sample mis-alignment. In any 
case, uncertainty in sample alignment has been checked 
to fall within 0.5°, and we now report angle resolved data, 
which reflect well the above finding, the small uncertainty 
in sample alignment notwithstanding. 

Fig.'s [2][a,b) show representative EDCs at (0, 0.5) 
(throughout this paper, we use 7r/a as unit of momen- 
tum, where a is the lattice constant of the Cu02 plane). 
Note that the x = 0.16 sample shows a clear isotope 
shift, ss 30 meV, while the x — 0.18 sample shows no 
isotope shift within the uncertainty ±10 meV. Panels 
c,d show the effect of isotope substitution on the posi- 
tion of high energy peak in the momentum-energy plane. 
The cuts investigated here are along the (0,0) to (0,1) 
direction, where we expect large isotope effectP^. Note 
that the peak positions of EDCs for optimal doping show 
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a definite isotope-induced energy shift as all of the red 
symbols are above the blue symbols, with an exception of 
one data point in cut 2. The shift is large, up to 40 meV, 
in good agreement with our previous worlpl^l. The effect 
is observed also for super-structure replicas (cuts 2,3), for 
positive and negative k y values (cut 2), and for k values 
in the first and the second Brillouin zones (cut 3). That 
is, this anomalous isotope effect observed has the correct 
expected symmetry. Also, it can be noted that as the 
binding energy is increased beyond 300 meV, the isotope 
effect decreases. Thus, the "active energy range" of the 
anomalous isotope effect is 200 — 300 meV, consistent 
with our previous resultP. In contrast, the over-doped 
sample shows no isotope effect within the error bar, as 
shown in panel cP^. 

The data presented above disagree with the data 
of Douglas et alP, assuming their sample is optimally 
doped. In view of the sensitivity of the isotope effect to 
a minute doping change, we decided to take a closer look 
at the doping value of their samples. 

Fig. [3] shows that the doping of samples in Ref. [I] is 
significantly different from ours. This can be seen from 
the comparison of the dispersions along the nodal cut, 
as shown in the main panel. In an earlier comprehensive 
studjE] it has been shown that the ratio of the high en- 
ergy (/i — E = 100 — 300 meV) group velocity to the low 
energy (/i — E < 50 meV) group velocity is a sensitive 
indicator of the doping level, reflecting the fact that the 
ARPES "kink" becomes milder as the doping increases. 
This ratio for the data of Ref. 1 (as 1.5) is smaller than 
that for our data (ss 1.7) by ss 0.2, which corresponds^ 
to an over-doping by Ax — 0.02 ~ 0.03. 

An over-doping of the samples in Ref. [1] is also appar- 
ent from both the Fermi surface and the superconducting 
gap size. As the inset of Fig. [3] shows, their Fermi sur- 
face is very different from the well-known optimal doping 
Fermi surface and, actually, looks like that of an over- 
doped sample^ with x at least 0.22. Similarly, the ab- 
sence of a clear superconducting gap in their off-nodal 
dispersion dataP contradicts the well-known behavioral 
at optimal doping and is also indicative of over-doping. 

While the above observations make it certain that Dou- 
glas et al.'s samples are overdoped, some differences in 
data remain to be understood. For instance, note that 
the main cause of the difference in velocity ratio be- 
tween their data and our data in Fig. [3] is the rather 
large low energy group velocity in Douglas et al.'s data 
as comp ared to ou r data as well as other data in the 
literatur e 1 11 ! 12 ! 13 ! 1 -^, while the high energy group velocity 
is almost the same as that of our data. This is quite 
curious in view of the well-known factP^ that the high 
energy group velocity has a much greater dependence on 
doping, and may require a re-examination of the overall 
momentum scale. 

Lastly, note that the fact that our observed isotope 
effect in Fig. [3] vanishes as energy increases beyond 300 
meV also makes it difficult to ascribe such an effect to a 
sample misalignment or a doping difference between iso- 
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FIG. 3: The nodal cut data of Ref. U compared with our 
nodal cut data. Also we show line shape fit results f or a more 
extended energy range than we presented previousljP^. Inset 
compares our Fermi surface with the Fermi surface of Ref. [T] 

tope samples. In addition, the latter possibility is ruled 
out firmly by the analysis that we now present. 

Fig. [4] presents the detail of how we characterize the 
hole doping. It is important to note that in order to 
obtain the absolute x value from the Fermi surface, the 
ARPES data in the entire Brillouin zone is required. This 
is done as shown in Figjfl and our previous work (sup- 
plementary information/^ However, one should keep in 
mind that the absolute value of x is made uncertain by 
line shape issues, such as broad line shapes, (pseudo) 
gaps, bi-layer splitting, and photon energy dependence of 
spectral functions, all of which become increasingly im- 
portant as momentum approaches the anti-nodal region. 
Thus, it is not surprising that the value of x typically has 
an error bar of magnitude ±0.01 or larger^. 

On the other hand, it is possible to compare the relative 
values of doping between two samples, by concentrating 
on the nodal region, where these problems are minimal. 
Fig. Ka) compares the intensity profile along the nodal 
cut for our three different optimally doped samples used 
in Fig. [3] Also included is the data for a reference sam- 
ple with a high over-doping^. As the doping x increases 
the Fermi surface moves toward the T point. The key 
observation is that, as shown by white arrows in panel 
b, a doping change causes M, SI, and S2 to move in 
the same direction while it causes S2' to move in the 
opposite direction. Thus, it is possible to judge the rel- 
ative change of doping quite accurately by focusing on 
spectral peaks, e.g. peaks from SI and S2', in a narrow 
angular range near the normal emission - the least com- 
plicated procedure. For instance, peaks from SI and S2' 
get closer, while peaks S2 and S2' move farther apart, 
as doping increases (see white arrows in a, b). Panel 



4 




OPT92 0.1 
x = 0.16 



kit (1,1) (it/a) 



(b)' 












/ S1 \ 






-4S2\ . S 2 - 

























(C) — x=0.18 

=== !X =0.16 , 'S1 




S2 S2' ==ibo 



x=0.16 x=0.18 





k//(1,1)(it/a) 



0.3 



FIG. 4: (a) Momentum distribution curves at constant en- 
ergy (35 meV below fi in order to include the S2' feature 
nicely within the cut) for optimally doped samples and a ref- 
erence over-doped sample^, (b) Tight binding Fermi surface 
and its superstructure replicas for x = 0.16 and x = 0.18 (in- 
set). The same diagonal rectangle, enclosing central regions 
of SI, S2 and S2', is indicated in both plots. Within that 
rectangle, note subtle but clear relative position changes of 
SI, S2 and S2' on over-doping, i.e. SI and S2' move closer 
while S2 and S2' move farther apart, (c) ARPES intensity at 
H, momentum-integrated over the short axis of the rectangle 
in (b), for statistics, and plotted as a function of momentum 
along the long axis of that rectangle. The doping-induced 
changes expected from (b) are clearly observable. For this 
plot, SI and S2 features for different doping values are aligned 
by a small momentum shift. 



a illustrates how the two peaks (at binding energy 35 
meV) from SI and S2' almost merge with each other for 
x = 0.20, while the position of S2' is still discernible as 
a shoulder. Using this as reference, it is then possible 
to estimate how much the S2' peak would have moved 
when the doping increases by 0.01. This movement, in 
good agreement with the tight binding parametrizatioEP, 
is indicated in the lower left between the two curves for 
16 and 18 0. It is obvious that all three optimally doped 
samples shown have the same doping with an error bar 
well within 0.01. Fig.'s !4{b,c) show a similar analysis for 
the data of Fig.'s 1,2. Here, the evolution of the Fermi 
surface from x = 0.16 to x — 0.18 is again clearly re- 
flected in the relative position changes of the SI, S2, and 
S2' features. 

That a mere 2 % change of doping level can cause a 
significant change in the response to isotope exchange 
may be at first quite surprising. On the other hand, 
such spectroscopic sensitivity to doping has a precedent 
in another strongly correlated material, namely the man- 
ganites. There, a recent studj^ reports that a doping 
change by 0.02 causes a qualitative change in the elec- 
tronic structure, a complete disappearance of the quasi- 
particle peak, again indicating that a seemingly insignif- 
icant doping change can be very crucial. 

In conclusion, we presented data on optimally doped 
and slightly overdoped samples of oxygen isotope ex- 
changed Bi2Sr2CaCu208+<5 samples. Surprisingly the 
anomalous isotope effect observed in optimally doped 
samples disappear in overdoped samples. This study 
prompts a more careful study of the electronic structure 
near the optimal doping, where the mystery of the high 
temperature superconductivity seems to hide in more 
veils to be removed. 
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